Quantum Monte Carlo study of a dominant s-wave pairing symmetry in iron-based 

superconductors 
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We perform a systematic quantum Monte Carlo study of the pairing correlation in the S4 symmet- 
ric microscopic model for iron-based superconductors. It is found that the pairing with an extensive 
s-wave symmetry robustly dominates over other pairings at low temperature in reasonable param- 
eter region regardless of the change of Fermi surface topologies. The pairing susceptibility, the 
effective pairing interaction and the (it, 0) antiferromagnetic(AFM) correlation strongly increase as 
the on-site Coulomb interaction increases, indicating the importance of the effect of electron-electron 
correlation. Our non-biased numerical results provide a unified understanding of superconducting 
mechanism in iron-pnictides and iron-chalcogenides and demonstrate that the superconductivity is 
driven by strong electron-electron correlation effects. 
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A today's major challenge in the study of iron- 
based superconductors P-Ql is how to obtain an unified 
microscopic understanding of the different families of 
these materials, in particular, iron-pnicitides and iron- 
chalcogcnidcs, which distinguish themselves from each 
other with distinct Fermi surface topologies 0-0. In the 
past several years, the majority of the theoretical stud- 
ies of iron-based high temperature superconductors were 
based on models with complicated multi-d orbital band 
structures [8Hl8j| . The conclusions from these studies pro- 
vided a good understanding of iron-pnictides while draw- 
ing a very different picture regarding of the magnetism 
and superconductivity of iron-chalcogcnides because of 
their strong dependence on theoretical approximations 
and the topology of Fermi surfaces. Effective mod- 
els emphasizing local AFM exchange couplings appear 
to unify the understanding of superconducting states of 
both materials U 27 1 . However, they lack of a support 
from more fundamental microscopic electronic physics, 
and using some unbiased numerical techniques is believed 
to be the only opportunity to win this great challenge as 
Hartree-Fork type approaches are biased if the electronic 
correlation dominates in the system. 

Recently, it has been shown that the underlining elec- 
tronic structure in iron-based superconductors, which is 
responsible for superconductivity at low energy, is essen- 
tially governed by a two-orbital model with a $4 symme- 
try, the symmetry of the building block-the trilaycr FeAs 
or FeSe structure in iron-based superconductors. In the 
model, the dynamics of two S4 iso-spin components are 
weakly coupled so that the essential physics is controlled 
by a single S4 iso-spin component. Thus, a minimum 
effective model that captures the low energy electronic 
and magnetic properties can be well described by an ex- 
tended one-orbital Hubbard model near half- filling [2^]. 
Such a microscopic understanding provides a new oppor- 
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FIG. 1: (Color online)(a) Here red and white circles indicate 
different sub-lattice A and B. The dark solid lines indicate ti, 
dot blue lines indicates t 2 , and dash purple lines indicates t' 2 . 
(b) The energy band along high symmetry line in unfolded 
Brillouin Zone. Solid dark line: ti=0.3, t 2 — lA, t' 2 =-0.6; dash 
red line: ti=0.3, t 2 =1.2, i 2 =-0.8 and dot blue line: ti=0.8, 
£ 2 =1.2, t 2 =-0.8[2§]. (c) Fermi surface at half filling for ti = 
0.3, £2 = iA,t' 2 = —0.6 (a typical case for iron-pnictides [30I— 
Hi]), and (d) at electron filling < n >= 1.1 for ti = 0.8, t 2 = 
1.2, t 2 = — 0.8( a typical case for iron-chalcogenides )■ 



tunity to make use of highly controllable and unbiased 
numerical methods to study iron-based superconductors, 
in particular, to obtain a possible unified understanding 
of iron-pnicitides and iron-chalcogenides. 

In this Letter, we perform a systematic quantum 
Monte Carlo study of the pairing correlation in the S4 
symmetric microscopic model. We find that the pair- 
ing with an extensive Ai g s-wave symmetry robustly 
dominates over other pairings at low temperature in 
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FIG. 2: (Color online) Phase of the s x 2 +y 2, d x 2_ y 2, s xy and 



reasonable parameter region regardless of the change 
of Fermi surface topologies, the presence or absence of 
hole pockets at r point. For both iron-pnictides and 
iron-chalcogenides, the pairing susceptibility, the effec- 
tive pairing interaction and the (tt, 0) AFM correlation 
strongly increase as the on-site Coulomb interaction in- 
creases. Our study demonstrates that the superconduc- 
tivity in iron-based superconductors is driven by electron- 
electron correlation and the nesting between electron and 
hole pockets is not an essential physics in iron-based 
superconductors. This conclusion is behind many pro- 
posed effective models T^- 21 1 but could not be conclu- 
sively reached with approximated methods [3]. The fact 
that the extended s-wave is favored even in the case with- 
out hole pockets differs from the earlier simple conjecture 



ll|,ll2j. Our unbiased numerical results thus present 



a rather different picture form mean field approaches. 

As shown in Fig. [T] (a), the minimum extended Hub- 
bard model for a single S4 iso-spin component in the 
iron-square lattice is described by 

H = ti E( a icA+'^ + h.c.) 

rqcr 

icr icr 
icr icr 

+U E( n «t n a4 + n bi ^n b%i ) + fi^(n aia + n bla \l) 

i icr 

Here, a,i a (a| CT ) annihilates (creates) electrons at site 
with spin a (er=T, 4-) on sublattice A, bi a (b\ a ) annihilates 
(creates) electrons at the site with spin a (c=t, \) on 
sublattice B, n ala = a\ a a. la , n bia = b\ a b la , 77 = (±x,0) 
and (0, ±y). In the above model, for simplicity and clar- 
ity, we only keep a minimum set of parameters which 
include three key shortest hopping parameters that are 
responsible for the physical picture revealed by the S4 
symmetry [2^|. The selection of parameters in following 
studies does capture the essential physics of typical cases 
for iron- pnictides [lO-lil and iron-chalcogenides 0-0, as 
shown in Fig.l (b-d). 

Our numerical calculations were mainly performed on 
an 8 2 or a 12 2 lattice with periodic boundary conditions. 



The system was simulated using determinant quantum 
Monte Carlo (DQMC) at finite temperature. The basic 
strategy of DQMC is to express the partition function as 
a high-dimensional integral over a set of random auxiliary 
fields. The integral is then accomplished by Monte Carlo 
techniques. In our simulations, 8000 sweeps were used 
to equilibrate the system. An additional 45000 sweeps 
were then made, each of which generated a measurement. 
These measurements were split into fifteen bins which 
provide the basis of coarse-grain averages and errors were 
estimated based on standard deviations from the average. 
For more technique details we refer to Refs. 
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As magnetic excitation might play an important role in 
the superconducting (SC) mechanism of electronic corre- 
lated systems, we first studied the magnetic correlations 
in such system. We define the spin susceptibility in the 
z direction at zero frequency, 



x(q) = 



[ 9 dT E ^2e^-^{m id (T)-m jd ,(0))(2) 

■* d,d'=a,b i,j 



-Ht 



with 



-J 



m ia =ai f a n - a^a iir 



where m,i a {T)=e Hl m, a (0)e 

and mi b =b\fiit - b^b^. 

To investigate the SC property of iron-based supercon- 
ductors, we computed the pairing susceptibility, 



dr(AT(i,r)A a (3,0)) 1 



(3) 



where a stands for the pairing symmetry. Due to the con- 
straint of on-site Hubbard interaction in Eq. (1), pairing 
between two sublattices is favored and the corresponding 
order parameter A^(z) is defined as 

A^(j) = fa( Sl ^ a ^ bi + s ^ ~ a i|A+*«t) t = 
l 

with / a (<5i) being the form factor of pairing function. 
Here, the vectors <5i (2=1,2,3,4) denote the nearest neigh- 
bour (NN) inter sublattice connections where 5 is (±x, 0) 
and (0, ±y), or the next nearest neighbour (NNN) inner 
sublattice connections where 8' is ±(i,y) and ±(i, —y). 

We study four kinds of pairing form, as that sketched in 
Fig. [5] For s x 2 +y 2-~w&ve pairing, f s (8i) = 1. For d x 2_ y 2 
pairing, fd(8\) is 1 when <5i = (±x,0) and -1 otherwise. 

Another two interesting pairing forms are d xy -~wa,ve 
and extensive s xy -wave, 

d^-wavc : fd„ y {5[) = l(8[ = ±(x,y)) 

and f dxy (8' 1 ) = -l(8' l = ±(x,-y)), 
s^-wave :/,.„(*{) = 1, i = 1,2,3,4. (4) 

In Fig. El we present the spin susceptibility x(q) at 
different electron fillings and U for temperature T=l/6. 
A sharp peak at (it, 0) in Fig. [3] indicates the existence 
of AFM spin correlation in iron-based superconductors 
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FIG. 3: (Color online) Magnetic susceptibility x(l) versus 
momentum q on an 8 2 lattice. 
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FIG. 4: (Color online) Pairing susceptibility P a as a func- 
tion of temperature for different pairing symmetries at ti = 
0.3, t 2 = 1.4,4 = -0-6 (solid line) and tx = 0.3, t 2 = 1.2,4 = 
—0.8 (dash line) with U = 4.0 and < n > = 1.0 on an 8 2 lattice. 
In the inset: The effective pairing interaction P S[cy — P Sx „. 



close to half filling. From Fig. [31 one can notice that 
at half-filling, the peak is very sharp, and as the elec- 
tron filling < n > decreases from half filling, x(q) is 
reduced around the {tt, 0) point, which indicates that 
the AFM spin correlation is suppressed when the system 
is doped away from half filling. The (tt, 0) AFM spin 
fluctuations have been universally observed in all iron- 
based superconductors It is clear that our model 
and the results naturally provides an explanation for the 
stable AFM exch ang e couplings J 2 observed by neutron 



scattering [39|, |4l|, |42j in parental compounds. 

From the behavior of magnetic correlation shown in 
Fig.|3l it is also important to note that the (tt, 0) spin cor- 
relation does not depend on the presence of the hole pock- 
ets at T point, which indicates that such an AFM correla- 
tion is driven by electron-electron correlation rather than 
nesting between hole and electron Fermi pockets [lU . The 
(tt, 0) AFM correlation is also stabilized by the fact that 
the diagonal or NNN hopping parameters t 2 are larger 
than the NN hopping t\. As shown in Fig. [3l if t\ is 
significantly large, the (tt, tt) magnetic correlation can be 
dominant. As we will show later, in this case, d x i_ y i 
pairing can be significantly enhanced. 



FIG. 5: (Color online) The effective pairing interaction P Bxy 
P Sxy as a function of temperature for different U at t\ 
0.3, t 2 = 1.4,4 = -0.6 (solid line) and ti = 0.3, t 2 = 1.2,4 
—0.8 (dash line) on an 8 2 lattice. 



Fig. [4] shows the temperature dependence of pairing 
susceptibilities for different pairing symmetries. Within 
the parameter range investigated, the pairing susceptibil- 
ities for various pairing symmetries increase as the tem- 
perature is lowered. Most remarkably, s xy increases much 
faster than any other pairings symmetry at low temper- 
atures. This demonstrates that the s xy pairing symme- 
try is dominant over the other pairing symmetry near 
half filling. When other parameters are fixed, reducing 
4s = (^2+4)/2 from 0.4 to 0.2, or increasing the absolute 
value of 4 ■ one may also see that pairing susceptibilities 
with different symmetries are all enhanced, in particular, 
the s xy pairing susceptibility. 

In order to extract the effective pairing interaction 
in different pairing channels, the bubble contribution 
Pu(i,j) is also evaluated, which is achieved by replac- 
ing (a^bjfal+g^bj+s^) with (a^^^L^ A'+£t'T> in 
Eq. (3). In the inset of Fig. HI we plot P S;cy - P Siy for 
t\ = 0.3,^2 = 1-4,4 = —0.6 (solid line). It is apparent 
that P Sj . — P Sa . shows a very similar temperature depen- 
dence to that of P Sxy ■ The effective pairing interaction 
for P s , is found to take a positive value and to increase 
with lowering temperature. The positive effective pairing 
interaction indicates that there actually exists attraction 
for the s xy pairing. The effective s xy pairing interaction 
for ti = 0.3,4 = 1-2,4 = -°- 8 is also shown in FigH 
as dash line. Comparing results with different t' 2 , the ef- 
fective s xy pairing interaction is enhanced greatly as the 
absolute value of t' 2 increases. 

In Fig. [5l we present the effective pairing interaction 
as a function of temperature for P S[c at different U 



with ti = 0.3, h = 1.4,4 



-0.6 (solid line) and (b) 



ti = 0.3, <2 = 1-2,4 = —0.8 (dash line). One can see 
that, the effect pairing interaction is enhanced as U in- 
creases. Especially, the effective s xy pairing interaction 
shown in Fig. [5] tends to diverge in low temperatures, 
and the increasing U tends to promote such diverge. This 
demonstrates that the electron-electron correlation plays 
a key role in driving the superconductivity. 
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FIG. 6: (Color online) The effective paring interaction P Sxy — 
P s as a function of temperature for different electron fillings 
at fi = 0.3, t 2 = 1.4, t' 2 = -0.6 (solid lines) and ti = 0.3, t 2 = 
1.2, t' 2 = —0.8 (dash lines) on an 8 2 lattice. 



We also studied the temperature dependence of effec- 
tive pairing interaction at different electron fillings. In 
Fig. |51 the temperature dependence of effective pairing 
interaction is shown for < n >=1.0, 0.9 and 0.8 with 
ti = 0.3, t 2 = 1.4,4 = -O- 6 ( a ); and < n >=l-0, 1.1 and 
1.2 with t x = 0.3, t 2 = 1.2,4 = -0.8(b). Both Fig. El 
(a) and (b) show that, the effective paring interaction de- 
creases as the system is doped away from half filling. As 
shown in Figs. (3,4,6), the decrease of the peak at (7r,0) 
of spin susceptibility is correlated with the suppression 
of the pairing susceptibility. This directly confirms that 
the (7T, 0) AFM fluctuations favors the s xy pairing. 

Following, we reported the effect of t\ on pairing sym- 
metry. In Fig. [Jj the pairing behavior for different t\ is 
shown for < n >= 1.1. As t\ increases significantly, it is 
possible that the d x 2_ y 2 wave pairing becomes dominant. 
A global phase diagram of the leading pairing symmetry 
in the t\ / t 2 -\t' 2 \ / 1 2 plane, obtained by DQMC at temper- 
ature T = 1/6, is shown in the inset of Fig. [Jj This 
phase diagram indicates that the s xy is robust when t\ 
is small. The phase transition line in the inset is roughly 
corresponding to the hopping parameter setting when the 
dispersion of the band at Y which is originally responsi- 
ble for the hole pocket gets reversed. Since such a reverse 
does not take place in the band structure of iron-based 
superconductors [30j] . the s xy -w&vc thus is robust. The 
robustness is further shown in Fig. [7J where we reported 
DQMC calculations throughout the t\jt 2 vs t' 2 /t 2 phase 
diagram. In general, the hopping parameters in the £4 
model is t 2 > t\ > \t' 2 \, and the fitting to band structure 
shows that t\ is roughly 0At 2 and t' 2 is around 0.3^9 [28|. 

Finally, we present the pairing susceptibility for dif- 
ferent pairing symmetries on a 12 2 lattice and compare 
them to the s xy paring susceptibility on an 8 2 lattice in 
Fig. [5] One can see that, the results for a larger lattice 
size confirm that the s xy dominates over other kinds of 
pairing symmetry. Moreover, it is interesting to see that, 
the s xy pairing susceptibility increases as the lattice size 
increases, especially at low temperature. This enhance- 
ment is consistent with the behavior of spin susceptibility 
x(q) shown in the inset of Fig. [51 in which x(q) is also 
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FIG. 7: (Color online) P a as a function of temperature at 
< n >=1.1, U = 4.0, t 2 = 1.2, t' 2 = -0.8 for different ti on 
an 8 2 lattice. Inset: the competition between s xy and d x 2_ y 2 
depends on ti/t 2 and —t' 2 /t 2 . 
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FIG. 8: (Color online) P a as a function of temperature for 
different pairing symmetries at ti = 0.3, t 2 = 1.4, t' 2 = —0.6 
of a 12 2 (dash line) and an 8 2 (solid line) lattice with U = 4.0 
and < n >=1.0. In the inset: x(l) versus q for a 12 2 (dash 
line with solid symbols ) and 8 2 (solid line with open symbols) 
lattice at < 71 >= 1.0 (dark red) and < n >= 0.9 (pink). 



enhanced as the lattice size increases. The peak of x( ( l) 
at (71", 0) point increases as the lattice size increases at 
< n >= 1.0, which implies there may be a static mag- 
netic order develops for the ground state near half filling. 
The above results of the pairing susceptibility have also 
shown to be valid when long range pairing correlation 
are calculated (45j. Thus, the 8 2 lattice is large enough to 
investigate the dominant pairing symmetry 

Overall, these results clearly suggest that the super- 
conductivity and pairing symmetry in iron-based super- 
conductors are determined by the combination of strong 
electron-electron correlation and the microscopic setting 
of hopping parameters. In the model, if t 2 and t' 2 are 
fixed, a small decrease of t\ can cause the vanishing of 
hole pockets at Y. Since the ci-wave pairing channel is 
caused by t\, we can conclude that the s-wave pairing 
is ever more robust in iron-chalcogenides than in iron- 
pnictides, a result completely different from weak cou- 
pling approaches. Recent experimental results by angle 
resolved photoemission spectroscopy have strongly sug- 
gested that pairing symmetry in both iron-pnictides [30| 
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and iron-chalcogcnidcs is a s-wave[46j. Our study clearly 
provides such a unified microscopic understanding. 

In summary, we study the paring susceptibility and 
effective paring interaction in iron-based superconduc- 
tors based on an effective S4 model [2a|, It is confirmed 
that the (ir, 0) AFM dominates at half filling and the 
paring susceptibility with s xy symmetry are enhanced 
as the electron-electron correlation increases, especially 
at low temperature. It is suggested that the s xy wave 
paring symmetry, which is the A\ g phase correspond- 
ing to the point group of the lattice and homogenous 
with respect to the transitional symmetry of the iron 
lattice4 47j , should be dominant in iron- based supercon- 
ductors. The reported strong AFM at half filling and 
the behavior of pairing susceptibility and effective pair- 
ing interaction strongly support that the microscopic su- 
perconducting mechanism for cuprates and iron-based 
superconductors including both iron-pnictides and iron- 
chalcogenides are identical. 
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